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Abstract 

We report on the multi-epoch observations of H2O maser emission in star forming region OH 43.8— OA 
carried out with VLBI Exploration of Radio Astrometry (VERA). The large-scale maser distributions 
obtained by single-beam VLBI mapping reveal new maser spots scattered in area of 0.7 x 1.0 arcsec, in 
addition to a 'shell-like' structure with a scale of 0.3 x 0.5 arcsec which was previously mapped by Downes 
et al.(1979). Proper motions are also obtained for 43 spots based on 5-cpoch monitoring with a time span 
of 281 days. The distributions of proper motions show a systematic outflow in the north-south direction 
with an expansion velocity of ~ 8 km s _1 , and overall distributions of maser spots as well as proper motions 
are better represented by a bipolar flow plus a central maser cluster with a complex structure, rather than 
a shell with a uniform expansion such as those found in Cep A R5 and W75N VLA2. The distance to OH 
43.8—0.1 is also estimated based on the statistical parallax, yielding D = 2.8 ± 0.5 kpc. This distance is 
consistent with a near kinematic distance and rules out a far kinematics distance (~ 9 kpc), and the LSR 
velocity of OH 43.8—0.1 combined with the distance provides a constraint on the flatness of the galactic 
rotation curve, that there is no systematic difference in rotation speeds at the Sun and at the position of 
OH 43.8—0.1, which is located at the galacto-centric radius of ~ 6.3 kpc. 

Key words: ISM:star forming regions — ISM:individual(OH 43.8-0.1) — masers (H 2 0) — VERA - 
Galaxy: kinematics and dynamics 



1. Introduction 

Star forming regions are often associated with wa- 
ter maser emission, and multi-epoch VLBI observations 
of water masers provide us a unique tool to study the 
structure and kinematics of star forming regions through 
proper motions of maser spots. Until the middle of 1990s, 
VLBI studies of H2O maser sources were made to reveal 
the kinematics in the well-known star forming regions such 
as Orion KL, W51M&N, Sgr B2, W49N, W3(OH) (Genzel 
et al. 1981a; 1981b; Schneps et al. 1981; Reid et al. 1988; 
Gwinn et al. 1992; Alcolea et al. 1993). These studies 
successfully detected proper motions of maser spots, and 
demonstrated that maser-emitting gases trace the shock 
regions where the outflows from forming stars hit ambient 
gases. Proper motions of maser spots were also used to de- 



termine the source distance via statistical parallax and/or 
detailed kinematic modeling, providing a powerful tool to 
determine distances in galactic scale, even up to ~ 10 kpc. 
Recently, numerous efforts were made to perform VLBI 
studies of more number of star forming regions, detecting 
proper motions of H 2 masers in several other sources 
such as IRAS 05413-0104, IRAS 21391+5802, S106 FIR, 
W3 IRS5, W51S, Cep A, W75N, NGC 2071, AFGL 5142, 
IRAS 20050+2720 (e.g., Claussen et al. 1998; Patel et al. 
2000; Furuya et al. 2000; Imai et al. 2000; 2002; Torrefies 
et al. 2001; 2003; Seth et al. 2002; Goddi et al. 2004; 
Furuya et al. 2005). While most of H 2 maser sources 
show bipolarity in their distributions as well as motions, 
there are a few exceptions exhibiting coherently expand- 
ing shell (Cep A R5 & W75N VLA2: Torrefies et al. 2001; 
2003), casting a question on the formation mechanism of 
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these morphologically different outflows in star forming 
regions. 

OH 43.8-0.1 (also referred to as G 43.8-0.1 or G 43.795- 
0.127) is a star forming region located in the galactic plane 
with OH, CH 3 OH and H 2 maser emissions. While op- 
tical and near IR observations tell little about this region 
due to the extremely large galactic extinction, this star 
forming region is known to be associated with an infrared 
source IRAS 19095+0930 and an Ultra-Compact HII re- 
gion (hereafter UC HII; Kurtz et al. 1994). VLA stud- 
ies of OH maser distributions revealed that OH maser 
spots are surrounding the UC HII (Argon et al. 2000). 
Molecular lines such as CO and CS were also detected 
toward the direction of OH 43.8-0.1 (Plume et al. 1992; 
Shepherd & Churchwell 1996), providing a systemic veloc- 
ity of Vlsr ~ 44 km s _1 . Assuming the standard galac- 
tic rotation (a flat rotation curve with R = 8.5 kpc and 
Oo = 220 km s _1 ), the systemic velocity gives a near kine- 
matic distance of 3 kpc and a far distance of 9 kpc. H2O 
maser emission in OH 43.8—0.1 was discovered by Gcnzcl 
& Downes (1977) and Dickinson et al.(1978), and shortly 
later VLBI observations of H2O maser were carried out 
by Downes et al.(1979). The distribution of maser spots 
obtained by Downes et al. (1979) suggested a shell-like 
structure with a diameter of 0.5 arcsec. Therefore, OH 
43.8—0.1 is another candidate for H2O maser sources that 
show shell-like morphology such as those recently found in 
Cep A R5 and W75N VLA2 (Torrelles et al. 2001; 2003). 
In order to confirm (or rule out) this, new multi-epoch 
VLBI observations are necessary to measure the proper 
motions of maser spots. 

Here we present results of multi-epoch observations OH 
43.8—0.1 H2O masers with VLBI Exploration of Radio 
Astrometry (VERA) . which is a new Japanese VLBI array 
dedicated to phase-referencing VLBI astrometry (Honma 
et al. 2000; Kobayashi et al. 2003). Originally the 
observations of OH 43.8—0.1 were carried out as a se- 
ries of test observations to evaluate the system perfor- 
mance of VERA, by observing OH 43.8—0.1 and an- 
other bright maser source, W49N, simultaneously in the 
dual-beam mode. The series of test observations already 
demonstrated a high capability of phase-referencing with 
VERA's dual-beam system (Honma et al. 2003), and 
also produced maser maps of W49N, successfully identi- 
fying the outburst spots which flared up in 2003 October 
(Honma et al. 2004). In the present paper, we focus 
on the results on OH 43.8—0.1, analyzed as single-beam 
VLBI observations. We present the maser spots map in 
larger scale than that of Downes et al.(1979), and measure 
proper motions for the first time for OH 43.8—0.1 based 
on 5-epoch observations spanning 281 days. We also de- 
termine the distance to OH 43.8—0.1 based on statistical 
parallax, and discuss the structure and kinematics of OH 
43.8—0.1 maser as well as constraint on galactic rotation 
curve obtained from the distance to OH 43.8—0.1. 



2. Observations and Reductions 

VERA observations of OH 43.8—0.1 were performed 
for 9 epochs between October in 2003 and September 
in 2004 with typical intervals of 1 ~ 2 months as a se- 
ries of test observations of VERA. The observations were 
done in dual-beam mode observing two bright H2O maser 
sources W49N and OH 43.8-0.1 (0°.65 separation) simul- 
taneously to evaluate the capability of phase-referencing 
as well as astromctric precision. The results of phase- 
referencing performance tests as well as dual-beam astro- 
mctric capability will be discussed elsewhere, and here we 
concentrate on OH 43.8—0.1. Left-hand circular polariza- 
tion signals were digitized and recorded using the VSOP- 
terminal system at a data rate of 128 Mbps with 2-bit 
quantization. Among the total bandwidth of 32 MHz, 
one 16 MHz channel was assigned to OH 43.8—0.1, cover- 
ing the radial velocity span of 215.7 km s _1 (adopting the 
rest frequency of 22.235080 GHz for H 2 6i 6 -5 23 tran- 
sition). Correlation processing was performed with the 
Mitaka FX correlator with a spectral resolution of 512 
points per channel, yielding frequency and velocity reso- 
lutions of 31.25 kHz and 0.42 km s _1 , respectively. 

In some epochs, one or a few stations were partly or 
fully missed due to system trouble and/or bad weather. 
Since VERA array has only 4 stations (Mizusawa, Iriki, 
Ogasawara, and Ishigaki-jima), lack of a station causes se- 
vere degrading of synthesized images, and would introduce 
a large error in doing multi-epoch measurements of maser 
proper motions. Hence, in the present paper, we use the 
data of 5 epochs that were obtained with the full 4-station 
array under relatively good conditions. The epochs pre- 
sented here are: day of year (DOY) 356 in 2003, DOY 
054, 120, 205, and 272 in 2004 (Dec. 22 in 2003, Feb. 23, 
Apr. 29, Jun. 23, and Sep. 28 in 2004). The system noise 
temperatures at the zenith were typically 150 K to 250 K 
for the first three epochs, and 200 K to 300 K for the last 
two epochs. 

In the data analysis of each epoch, calibrations of clock 
parameters, bandpass, amplitude and phase were carried 
out in a standard manner. Clock parameters (clock off- 
set and clock rate offset) were calibrated using the resid- 
ual delay and delay rate for a bright calibrator source, 
TXS 1923+210 (= ICRF J192559.6+210626), which was 
observed every 2 hours. The amplitude calibrations were 
done using the system noise temperature, which were eval- 
uated by the 'R-Sky' method, observing the reference 
black body at the beginning of each scan (typically ev- 
ery hour). Bandpass calibrations were done using auto- 
correlation spectra of continuum sources. For phase cal- 
ibrations, visibilities of all velocity channels were phase- 
referenced to the reference maser spot at Vlsr of 38.2 km 
s _1 , which is one of the brightest spots, and shows no 
sign of structure according to the closure phase. After 
those processes, calibrated visibilities were Fourier trans- 
formed to synthesized images, and the positions of the 
brightness peaks were determined with respect to the ref- 
erence maser spot (note that the reference spot was the 
same for all epochs) . The synthesized beam has a FWHM 
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beam size of 1.2x2.3 mas with a PA of 45°. In order to 
detect maser spots that are far from the reference spot, 
spot positions were first estimated using fringe rates, and 
the detailed mapping is performed around the estimated 
positions. 

For analyses of multi-epoch VLBI data, we identified 
maser spots based on spot positions as well as radial ve- 
locities. Spots in different epochs were identified as 'same' 
if their positions were within 1 mas and their radial ve- 
locities were within 0.21 km s _1 , which is a half of ve- 
locity resolution. For maser spots that were detected for 
more than 3 epochs, their proper motions were determined 
based on the linear least-squares-fit. 

3. Results 

3.1. Spectral Evolution 

Figure 1 shows the spectral evolutions of OH 43.8—0.1 
for the 5 epochs. The LSR velocities of maser lines range 
from 26 to 53 km s _1 , showing no major difference from 
the observations in 1977 (Downcs ct al. 1979). As seen in 
figure 1, auto-correlation spectra of OH 43.8—0.1 showed 
rapid variability of maser intensity during our observa- 
tions. The strong variability of OH 43.8—0.1 was already 
found by long-term monitoring by Lekht (2000), and sim- 
ilar variabilities are found in many of H2O maser sources, 
including rapid flare-up of a certain feature. The most 
prominent variation occurred at Vlsr of 38.2 km s _1 , 
which had l.OxlO 3 Jy on DOY 356 in 2003 and decreased 
to 75 Jy on DOY of 272 in 2004. Lekht (2000) also found 
a rapid flare at Vlsr =38.2 km s _1 in 1997 exceeding 
3.7xl0 3 Jy with a duration of 6 months. 

Figure 1 also shows the cross-correlated spectra of OH 
43.8—0.1 for Mizusawa — Ishigaki-jima baseline, which is 
the longest among the 6 baselines of VERA. The cross- 
correlation spectra in figure 1 were taken around the max- 
imum UV distance, being ~ 2200 km or 1.6xl0 8 A. In 
figure 1, for some features correlated flux is much smaller 
than total power flux, indicating that maser features are 
partially or fully resolved out with this baseline. The LSR 
velocities of cross-correlated maser lines mostly range from 
34 to 51 km s _1 , while some features around 25 to 34 km 
s _1 are detected in cross-correlated spectrum on DOY 120 
in 2004. The component of Vlsr =38.2 km s _1 is strongest 
in cross-correlated spectra in the first three epochs, with 
its peak intensity of 870 Jy on DOY 356 in 2003. The 
high ratio of correlated flux to total flux (8.7xl0 2 Jy to 
1.0x10 s Jy) suggests that the flaring of the 38.2 km s _1 
spot occurred in a compact area, being smaller than 6 AU 
(corresponding to 2 mas at the distance of 2.8 kpc; see 
subsection 3.4 for distance determination). In the follow- 
ing analyses, the maser spot at 38.2 km s _1 was taken to 
be the reference spot for visibility phase calibration. 

3.2. Maser Spot Distributions 

Figure 2 (a) shows the composite map of OH 43.8—0.1 
maser spot distributions with a size of 1 ".5 x 1 ".5. Spots 
that are detected at least in one epoch are plotted. The 
number of detected maser spots in each epoch are 53, 54, 
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Fig. 1. Five-epoch spectral evolution of OH 43.8—0.1 H2O 
masers. Thin lines show total flux (auto-correlation), 
and dashed lines show cross-correlated flux for 
Mizusawa — Ishigaki-jima baseline. The vertical arrow 
shows the flux scale of 500 Jy, and each spectrum is displaced 
by the same amount. 

60, 47 and 51 from DOY 356 in 2003 to DOY 272 in 
2004. For masers spots detected for more than 3 epochs, 
proper motion vectors are superposed (for proper motion 
measurements, see next subsection). 

The map origin of figure 2 (a) is the position of the 
reference maser spot at Vlsr =38.2 km s _1 . The ab- 
solute position of the map center is estimated to be 
(aj2ooo, <5 J 2ooo)=(19' I ll m 53 s .986, +09°35'50".22) with an 
uncertainty of .2 using the fringe rates of the reference 
maser spot. Color index in figure 2 shows the LSR veloc- 
ity range from 34.4 km s _1 to 50.9 km s _1 , where most of 
maser spots are located. Only exceptions are those with 
LSR velocities of 25 ~ 34 km s _1 , which were only de- 
tected on DOY 120 of 2004 and shown as open circles in 
figure 2 (b) and 2 (c). 

Figure 2 (b) is an expanded map of the shell-like struc- 
ture discovered by Downes et al. (1979). While figure 2 
(b) is consistent with the map obtained by Downcs ct al., 
figure 2 (a) shows that the shell-like structure found pre- 
viously do not represent the entire region of maser spot 
distributions. In fact, in figure 2 (a), the maser distribu- 
tions show an extension with a scale of ~ 1 arcsec from 



4 



M. Honma et al. 



the northern top of 'shell-like' structure toward the south. 
This extension of maser spots were not reported in previ- 
ous studies and discovered for the first time in the present 
study. The northern maser spots at (-145, 136) have Vlsr 
of 47 to 51 km s -1 , and the southern-most maser spot at 
(129, -810) has Vlsr of 41 km s" 1 . These properties of 
spot positions as well as radial velocities indicate that the 
maser spot distribution in OH 43.8—0.1 is not represented 
by a simple shell, but may be interpreted as a bipolar 
structure in the north-south direction. 



In addition to the north-south extension, there are a few 
other maser spots which are rather isolated, such as the 
one at (X,F) = (399, -107) and (-268, -507), which were 
also not reported in the previous studies. 




Fig. 2. Composite map of maser spot distributions and 
proper motion vectors of OH 43.8—0.1. Spot color indicates 
the LSR velocity (see color index at the top). The systemic 
velocity of OH 43.8-0.1 is Vlsr ~ 42 km s _1 . Left (a): the 
large scale map with a size of 1.5 arcsec X 1.5 arcsec (1 arc- 
sec corresponds to 2800 AU at the distance of 2.8 kpc). The 
arrow at bottom-right corner indicates the proper motion of 
0.5 mas yr~ x , corresponding to 6.6 km s _1 at the distance of 
2.8 kpc. Top-right (b): expansions of the shell-like structure 
found by Downcs et al.(1979). Bottom-right (c): close-up to 
the central cluster with proper motion vectors. Circles in (b) 
and (c) are maser spots with radial velocities out of the color 
index range (Vlsr=25 to 34 km s" 1 ) detected only in DOY 
120 of 2004. 



3.3. Proper Motion Measurements 

As described in the previous section, maser spots in 
different epochs were identified as 'same' if their po- 
sitions and velocities are within 1 mas and 0.21 km 
s -1 , respectively. Based on these identification criteria, 
43 maser spots were identified in more than 3 epochs, 
and proper motions were determined based on the lin- 
ear least-squares-fit. Table 1 summarizes the results of 
proper motion measurements obtained by a linear least- 
square-fit, including LSR velocities, spot positions, proper 
motions, detected epochs, and RMS residuals. Note 
that the RMS residuals 6x and Sy are defined as 5x = 
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yS^Xobg — Xat) 2 /n, where X Q ^ S and Xg t are observed 
and best-fit positions of spots and n is the number of data 
(3 ^ n S 5). Values of Sx and Sy provide us uncertainty 
estimates in the proper motion measurements. The max- 
imum values of Sx and Sy are 0.060 and 0.082 mas, and 
the mean values of Sx and Sy for 43 spots are 0.022 and 
0.029 mas. These values indicate that the proper motion 
measurements were successfully done with an accuracy 
better than 0.1 mas. Also, in table 1 there is no system- 
atic trend in Sx and Sy with LSR velocity or positions 
of maser spots. Since the maser spots are bright enough 
(typically larger than 10 Jy in correlated flux), the resid- 
uals from the best-fit linear motion are likely to be dom- 
inated by internal structure of maser spots rather than 
thermal noise. In fact, adopting the distance of 2.8 kpc 
(see next subsection for distance determination), the mean 
value of Sx and Sy correspond to 0.06 AU and 0.08 AU. 
Observations of other H 2 maser sources suggested that 
maser clouds have internal structure with a typical scale 
of 1 AU and a typical velocity width of 0.5 km s _1 (e.g., 
Gwinn 1994). Assuming the same values for masers in 
OH 43.8—0.1, one can roughly estimate possible changes 
in maser peak positions due to internal turbulent motions 
during our time span of observations (281 days) as 0.5 km 
s _1 x 281 x 86400 s = 0.08 AU, which is quite close to the 
physical scales given by the mean values of Sx and Sy . 

The proper motion vectors of maser spots listed in ta- 
ble 1 are shown in figure 2 (a). Proper motion vectors 
are shown for the spots with proper motions greater than 
0.1 mas yr _1 (40 among the 43 spots in table 1). As de- 
scribed in the previous section, the maser spot at 38.2 km 
s _1 was taken to be the reference maser spot, and thus 
all the proper motion vectors arc relative to the proper 
motion of the reference spot itself. The proper motion 
of the reference maser spot could not be measured in the 
present study, and thus it might introduce a bias to proper 
motions of other maser spots if the reference maser spot 
itself has a large proper motion. However, the mean val- 
ues of nx and \xy for 43 spots in table 1 are 0.117 mas 
yr _1 and 0.086 mas yr _1 , corresponding to 1.6 km s _1 
and 1.1 km s _1 at the source distance of 2.8 kpc. These 
values of mean proper motions are fairly smaller than the 
radial velocity width (~25 km s _1 , see figure 1) and typ- 
ical proper motion in table 1 (~0.7 mas yr -1 ), indicating 
that the proper motion of the reference maser spot does 
not introduce large systematic motions of the reference 
frame in figure 2. 

As seen in figure 2 (a), the global distribution of proper 
motions show bipolar structure with red-shifted northern 
components around (X. V) = (-145, 136) moving toward 
the north and moderately blue-shifted southern compo- 
nent at (129, -810) moving toward the south. Note that 
while there are more than 10 spots in the red-shifted 
northern components, there is only one blue-shifted south- 
ern component for which proper motion was measured. 
This is likely to be due to the geometric effect of bipo- 
lar flow which causes the difference in the column density 
of masing gas toward these components: red-shifted com- 
ponents locate on the far side of the protostar along the 



line-of-sight and thus the column density of masing gas is 
higher than the blue-shifted components, which locates on 
the near side of the protostar. Between these two compo- 
nents there is a cluster of maser spots with X between -80 
and -20 and Y between -310 and -250, which hereafter we 
refer to as the 'central cluster', shown in figure 2 (c). The 
directions and amplitudes of proper motion vectors in the 
central cluster are quite different from those of northern 
and southern components, showing rather complex struc- 
tures. For instance, maser features at (X, F)=(-83, -276) 
(spot number 13, 15 to 18 in table 1) and at (-18, -310) 
(spot number 29, 31 to 33) exhibit large velocity changes 
within the maser features in 0.3 mas scale (~ 0.8 AU). 
Highly blue-shifted maser spots at radial velocity of 25 to 
34 km s _1 (which were detected only for DOY 120, 2004) 
also locate in the central cluster area (open circles in figure 
2 c) . These complex structure of the central cluster may 
indicate that the protostar is associated with the central 
cluster, and the protostar's activity causes the complexity. 

In contrast to an expanding shell such as those found in 
Cep A R5 (Torrefies et al. 2001) or W75N VLA2 (Torrefies 
et al. 2003) , the proper motions in the shell-like structure 
found by Downes et al. cannot be represented by a coher- 
ent expansion from the shell center. The proper motions 
of the 'shell' in the north-south direction are 0.5 to 0.7 mas 
yr -1 , and are much larger than those of the east- west di- 
rection, which arc around 0.1 to 0.2 mas yr" 1 . The global 
properties of proper motion distributions indicate that the 
structure of OH 43.8—0.1 H 2 maser can be represented 
by bipolar flow plus central cluster, and that the shell-like 
structure in figure 2 (b) is just a part of bipolar flow and 
the central cluster (for further discussion on the structure 
of OH 43.8—0.1, see next section). 

3.4- Distance Determinations by Statistical Parallax 

Distances to maser sources can be determined from the 
dispersions of radial velocities and proper motions based 
on statistical parallax. Statistical parallax assumes only 
that the velocity dispersions are isotropic, and the basic 
equation for distance determinations with statistical par- 
allax can be written as, 
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Here D is the source distance, cry is the radial velocity 
dispersion, and ct m is the one-dimensional proper motion 
dispersion, and a^x and a^y are the proper motion dis- 
persions in X and Y directions, respectively. In practice, 
however, equation (1) is too much simplified, and one has 
to take into account the effect of observational uncertain- 
ties using following equation instead of equation (1), 
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The e terms correct for the effect of observational uncer- 
tainty that broadens the widths of velocity and proper 
motion. In the following analyses we make use of equa- 
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tion (2) as the basic equation for distance determination 
based on statistic parallax (but note that equations (1) 
and (2) give almost identical results because the e terms 
are sufficiently small in the present study). 

In several cases H 2 maser sources show bipolarity and 
thus anisotropy, and so careful analyses are necessary to 
obtain reliable distances. In the present paper, we perform 
statistical parallax analyses based on four different defini- 
tions of a^x and er^y : I) using all maser spots with disper- 
sion a defined as root-mean-square value from the average 
(i.e., a^x = -\/Wfix — ^x,i) 2 /N), 2) using all maser spots 
with dispersion a defined as the peak-to-pcak value (i.e., 
&(iX = max(//xi) — min(^xi); 3) the same definition of a 
to I) but only with spots in the central cluster, and 4) the 
same definition of a to 2) but only with spots in the cen- 
tral cluster. Note that in all cases, the dispersions of LSR 
velocities arc also defined in the same manner with those 
of proper motions (i.e., oy = \fT,{V — Vi) 2 /N for cases I) 
and 3), and ay = max(T^) — min(Vi) for cases 2) and 4). 
Cases 3) and 4) are to avoid the anisotropy caused by the 
bipolarity of maser flow: as already seen in the previous 
sections, while the large scale proper motions are repre- 
sented by a bipolar flow, the proper motions in the central 
cluster show rather complex structure. Also, cases 2) and 
4) are for testing how different definitions of 'dispersion' 
would affect the results, as the distributions of LSR veloc- 
ities and proper motions are not necessarily modeled by a 
Gaussian distribution. In all cases, ey is taken to be 0.21 
km s" 1 , which is the half of velocity resolution and the 
threshold for identifying maser spots in different epochs, 
and ex and ey are taken to be 0.029 and 0.038 mas yr" 1 , 
which are determined from averaged values of 5x and 5y 
and the time span r=281 day=0.77 yr, as ex = $x/ T and 
ey = Sy/r. 

Table 2 summarizes the results of distance determina- 
tion with case 1) to 4). Calculated distances vary from 
3.3 kpc for case 1) to 2.3 kpc for case 4). There are also 
some systematic trends in 4 cases: analyses with all spots 
tend to give larger distances, and analyses with RMS def- 
inition of dispersions tend to give larger distances. The 
former tendency may reflect a geometric effect of bipolar- 
ity. For instance, in cases 1) and 3) dispersions of proper 
motions are similar, but dispersions of LSR velocities are 
increased by factor of 1.26 when all spots are included 
(case 1), yielding the farther distance. This result may 
imply that the location of the bipolar structure is head- 
on, i.e., the angle between bipolar axis and line-of-sight 
is less than 45 degrees. The origin of latter tendency in 
table 2 (larger distance for RMS definition of dispersions) 
is unclear, but the difference of distances for two defini- 
tions of a is a factor of 1.14, and this could be an error 
introduced by small number statistics (7V=16 or 43). 

It is not easy to judge which of the four cases best 
represents the true distance to OH 43.8—0.1, as all four 
cases have advantages and disadvantages. While analyses 
with spots only in the central cluster avoid the effect of 
bipolarity, the number of spots become smaller and de- 
grade statistic reliability. The peak-to-pcak definition of 
dispersions is severely affected by an outlier with an ex- 
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Table 2. 



[Vol. 



Results of distance determinations 



case 


spot a 


N b 


defer c 


a' d 




D 1 


1 


all 


43 


RMS 


4.69 


0.303 


3.28 


2 


all 


43 


PP 


16.5 


1.23 


2.82 


3 


CC 


16 


RMS 


3.72 


0.294 


2.67 


4 


cc 


16 


PP 


10.1 


0.930 


2.30 












mean 


2.8 ±0.5 



a: spot selection: all for all spots, and CC for spots in 

the central cluster. 

b: number of spots used in analyses 

c: definitions for cry, a^x, and cr„y: RMS for root-mean- 
square, and PP for peak-to-peak definitions (see text). 



d: a'y in km s _1 (equation [2]). 
e: cr' in mas yr" 1 (equation [2]). 
f: determined distance in kpc. 



treme LSR velocity or proper motions (if exists), while 
the RMS definitions of dispersions may not work well for 
non-Gaussian distributions (e.g., Vlsr distribution in fig- 
ure 1). In order to obtain most-likely value of the distance 
to OH 43.8—0.1, here we conservatively take an average 
of 4 cases, which yields D = 2.8 ± 0.5 kpc. This distance 
is consistent with the near kinematics distance (~ 3 kpc), 
and rules out the far kinematic distance (~ 9 kpc). 

4. Discussion 

4.1. Structure of OH 43.8- 0.1 H 2 masers 

The key to understand the structure and kinematics 
of H2O maser in the context of star formation in OH 
43.8—0.1 is to know the location of the forming star rela- 
tive to maser spots in figure 2. Unfortunately H2O maser 
spots are most likely to be located in shock regions where 
outflow from a forming star hit ambient gases, and thus 
the position of the forming star cannot be directly ob- 
tained from H2O maser observations. An approach to ob- 
tain the position of a forming star is to use the continuum 
radio emissions, as a forming OB star radiate strong UV 
photons and produce a compact HII region. In fact, an 
Ultra-Compact HII region (UC HII) was found by VLA 
observations at frequencies of 8 GHz and 15 GHz (Kurtz 
et al. 1994). They found an unresolved UCHII at the posi- 
tion of (ajaooo, <5,/2ooo)=(19' l n m 53 s .987, +09°35'50".308) 
with an uncertainty of 0.1. This agrees well with our 
absolute position determination for the reference spot, 
(19' l ll Tn 53 s .986, +09°35'50".22) ±0".2. Thus, the UC HII 
regions indeed locate in the maser emitting regions, but 
precise location in figure 2 (a) is still unclear due to posi- 
tion error. 

As for OH maser emissions, Argon et al.(2000) located 
OH maser spots with respect to UC HII region using VLA 
A- Array. Their results demonstrated the clear association 
of the OH maser with the UC HII region. It is notable that 
the LSR velocity of OH maser emissions ranges from 40 to 
45 km s _1 (Argon et al. 2000), which is fairly close to LSR 
velocity range of the central cluster shown in figure 2 (c) 
(34 to 45 km s" 1 ). The results of OH maser observations 
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by Argon et al. (2000) may imply a connection between 
OH maser and the central cluster of H2O maser, and thus 
a possible link between the central cluster and the UC HII 
region. 

In addition to OH maser, methanol (CH 3 OH) masers 
at 6.7 GHz were found in OH 43.8-0.1 (Menten 1991, 
Szymczak et al. 2000). The 6.7 GHz methanol maser 
is categorized as class II methanol maser, which shows a 
clear connection with star-formation, presumably associ- 
ated with accreting matter or accretion disk around a pro- 
tostar. Mcntcn (1991) reported that the methanol maser 
line in OH 43.8-0.1 has Vlsr range of 39 to 44 km s _1 . 
Szymczak et al. (2000) also confirmed this result, while 
their results show that the peak intensity occurs at 39 to 
40 km s _1 . These Vlsr ranges are closer to those of the 
central cluster (34 to 45 km s _1 ), rather than northern 
maser spots at (A, F)=(-145, 136) (46 to 51 km s" 1 ). 
Thus, a simple interpretation of Vlsr of methanol maser 
lines is that methanol maser emissions are associated with 
the central cluster. If these methanol emissions are also 
associated with a forming star and thus UC HII, this in- 
dicates that the unresolved UC HII is located in or near 
the central cluster. 

The other way to estimate the position of forming star 
with respect to maser features is to use kinematic informa- 
tion. Since the outflow is powered by the forming star, its 
position can be estimated as an origin of outflow motion. 
Detailed modeling have been done for several sources to 
evaluate the position of outflow origin as well as distances 
(e.g., Reid et al. 1988; Gwinn et al. 1992; Imai et al. 
2000) . Here we take one of the simplest (and often used) 
kinematic models which assumes a radial outflow from 
a single origin with a constant outflow velocity. Basic 
technique is essentially the same as the one used in Imai 
et al.(2000), where the outflow origin was determined to 
minimize the differences between observed velocity vec- 
tors and model velocity vectors. Based on this kinematic 
modeling, the outflow origin was estimated to be at (A, 
F)=(-260 ± 130, -250 ± 130) mas. Although the error 
bars are relatively large, the estimated position of the 
outflow origin coincides with the position of the central 
cluster within the error bars. Considering the discussions 
described above, one can infer that the central maser clus- 
ter is associated with a forming star and the bipolar flow 
in the north-south direction is powered by the forming 
star in or near the central cluster. This kind of structure, 
a bipolar flow plus a central cluster, is also seen in other 
H2O maser sources, such as IRAS 05413-0104 (Claussen et 
al. 1998) and IRAS 21391+5802 (Patel et al. 2000), which 
are low mass star forming regions. To obtain more firm 
conclusions on the case for OH 43.8—0.1, precise colloca- 
tion of H2O maser spots with UC HII regions is definitely 
required. 

The distance to OH 43.8—0.1 provide us a physical scale 
of H2O emitting regions. The north-south bipolar exten- 
sion (0.98 arcsec) corresponds to a physical scale of 2740 
AU at the distance of 2.8 kpc. This scale is remarkably 
smaller than the maser structures of other massive star 
forming regions such as W49N (~ 20000 AU, Gwinn et al. 



1992), Sgr B2(N) (~ 28000 AU, Reid et al. 1988), and 
Orion-KL (~ 19000 AU, Genzel et al. 1981a). The size 
difference mainly originate from the small outflow veloc- 
ity of OH 43.8—0.1 rather than age differences in these 
star forming regions. In fact, dynamical age estimate of 
OH 43.8-0.1 H 2 maser gives r dyn = 750 ± 100 yr for the 
north-south bipolar extension, which is not much different 
from those of massive star forming regions, for instance, 
r dyn ~ 500 for W49N (Gwinn et al.1992). On the other 
hand, H2O maser motion in these massive star forming re- 
gions exceed 100 km s — 1 , while that in OH 43.8—0.1 is ~ 
8 km s _1 . Thus, this scale difference reflects the difference 
of outflow energy between OH 43.8—0.1 and other mas- 
sive star forming regions, and this is presumably linked 
to lower power, implying lower mass of a forming star in 
OH 43.8-0.1 than the protostars in W49N, Sgr B2 and 
Orion-KL. 

4-2. Constraints on Galactic Rotation Curve 

Although this study presents a distance measurement 
of only one source, the result can be used to constrain 
the rotation properties of the Galaxy. Assuming a circu- 
lar rotation at any radius of the Galaxy, observed radial 
velocity can be written as, 



V r = 



I If J 



sin/. 



(3) 



Here R and Rq are the galacto-centric distances of the 
source (OH 43.8—0.1 for the present study) and the Sun, 
and Q(R) and Oo are circular rotation speed of the Galaxy 
at the radius of the source and the Sun, and I is the galac- 
tic longitude. Note that V r and I are observablcs, and R 
can also be regarded as an observable because R is written 
as 



R= \ D 2 + RZ-2DR cosl, 



(4) 



and hence one can immediately know R once D and Rq 
are given. Equation (3) contains three parameters to be 
determined, namely, Rq, Oo, 0(i?). Among them, the 
galactic constants, Rq and Oo, are not necessarily un- 
known as there have been numerous studies to evaluate 
them. Therefore, if Rq and Oo are given, then 0(i?) can 
be determined, giving a constraint on the galactic rotation 
speed. 

The IAU standard values are Rq = 8.5 kpc and ©0 = 220 
km s _1 (Kerr & Lynden-Bell 1986), but several recent 
studies claim different values. Determinations of Ro later 
than 1986 tend to give smaller values, and Reid (1993) 
summarized the determinations of Rq, giving Rq = 8.0 ± 
0.5 kpc. As for O the situation is much complicated: 
while number of recent studies with various methods gave 
values around 200 ~ 220 km , being in good agreement 
with previous ones (e.g., Mignard 2000; Bcdin et al. 2003; 
Kalirai et al. 2004; Reid and Brunthaler 2004), there 
were also some studies which report fairly smaller or larger 
values varying from 180 to 270 km s _1 (e.g., Oiling & 
Merrifield 1998; Miyamoto & Zhu 1998). In the following 
analyses, we conservatively take large error for both Rq 
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and 6 , taking R = 8.0 ± 1.0 kpc and 9 = 220 ± 50 km 
s- 1 . 

Here we rewrite equation (3) by introducing a new pa- 
rameter 8 r = 8(i?)/0o as, 

\B smi / Rq 

Based on the galactic constants Rq and Go, one can deter- 
mine 6 r , which is the ratio of the galactic rotation speeds 
at the source and at the Sun. If the rotation curve is com- 
pletely flat as found in many bright spiral galaxies (e.g., 
Sofuc & Rubin 2000), 9 r should be unity. In order to 
obtain 6 r from equation (5), one also have to know the 
value of V r , the radial velocity. Assuming that the central 
cluster of H 2 maser is indeed associated with the UC 
HII regions and thus its radial velocity reflect the system- 
atic velocity of OH 43.8—0.1, the averaged velocity give 
Vlsr = 40 km s _1 . Observations of molecular gases with 
CS and CO lines give Vlsr = 43 — 44 km s" 1 (Plume ct 
al. 1992; Shepherd k Churchwell 1996), and methanol 
maser lines and OH maser lines lie in the LSR velocity 
range of 39 to 44 km s _1 (Menten 1991; Szymczak et 
al. 2000; Argon et al. 2000). Summarizing these results, 
here we adopt V r = 42 ± 3 km s _1 as the systemic velocity 
of OH 43.8-0.1. Using these values (R Q = 8.0 ± 1.0 kpc, 
9 = 220±50 kms" 1 , L> = 2.8±0.5 kpc, and V r = 42 ±3 
km s _1 ), equation (5) yields that 

*r = 1.00t°;S°, (6) 

where the upper and lower limits are calculated as the 
worst values with possible combinations of above pa- 
rameter ranges. The value of 8 r indicates that there is 
no systematic difference between rotation speeds at OH 
43.8—0.1 and at the Sun, and thus that the galactic ro- 
tation curve is fairly flat, being consistent with previous 
studies (e.g., Honma & Sofue 1997). For R = 8.0 kpc, 
equation (4) gives R = 6.3 kpc as a galacto-centric dis- 
tance of OH 43.8—0.1, and the interpolation of rotation 
speed between the Sun and OH 43.8—0.1 suggests that 
flat rotation curve spans nearly 2 kpc. 

5. Summary and Concluding Remark 

We carried out 5-epoch observations of OH 43.8—0.1 
H2O maser with VERA, and successfully measured the 
proper motions of OH 43.8—0.1 H2O masers for the first 
time. The large-scale maser map revealed north-south ex- 
tension of maser spots, and proper motion measurements 
showed that the distribution as well as the kinematics of 
OH 43.8—0.1 H2O maser can be represented by the bipo- 
lar flow in north-south direction plus the central cluster 
which has rather complex velocity structure. Based on 
statistical parallax, we also determined the distance to 
OH 43.8—0.1 as 2.8 ± 0.5 kpc, ruling out a far kinematic 
distance. The distance, combined with previous studies 
on the galactic constants such as Rq and Oo, gave a con- 
straint on the flatness of galactic rotation curve, showing 
no systematic difference in the galactic rotation velocity 
at the Sun and OH 43.8-0.1. 



In the near future, VERA will observe hundreds of H2O 
maser sources to perform astrometry of galactic star form- 
ing regions. While the main target of VERA is to investi- 
gate the structure and dynamics of the Milky Way Galaxy 
via astrometry of maser sources, VERA will also provide 
useful information on internal structures and kinematics 
of star forming regions, and the present study is one of 
the first examples of such valuable by-products of VERA. 
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Analysis Center (ADAC) of the NAOJ, which is an inter- 
university research institute of astronomy operated by 
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Table 1. Results of proper motion measurements. 



ID Vlsr X Y_ MX Mr epoch S x $Y 



1 


50.92 


-145.105 


136.124 


0.129 


0.152 


00111 


0.012 


0.027 


2 


50.47 


-145.085 


136.030 


0.088 


0.416 


11111 


0.028 


0.035 


3 


50.04 


-145.132 


136.030 


0.202 


0.461 


11111 


0.031 


0.026 


4 


49.62 


-145.117 


135.997 


0.180 


0.529 


11111 


0.030 


0.028 


5 


49.20 


-145.102 


135.982 


0.180 


0.552 


11111 


0.017 


0.030 


6 


48.78 


-145.092 


135.961 


0.146 


0.689 


11111 


0.016 


0.024 


7 


48.36 


-145.089 


135.957 


0.146 


0.740 


11111 


0.021 


0.016 


8 


47.94 


-145.093 


135.998 


0.144 


0.699 


11111 


0.020 


0.036 


9 


47.52 


-145.084 


136.015 


0.140 


0.647 


11111 


0.015 


0.022 


10 


47.10 


-145.076 


136.030 


0.166 


0.570 


11111 


0.011 


0.023 


11 


46.67 


-145.090 


136.040 


0.197 


0.518 


11111 


0.026 


0.039 


12 


46.29 


-155.368 


134.127 


0.088 


0.325 


10101 


0.035 


0.029 


13 


45.02 


-82.968 


-276.683 


-0.011 


0.204 


00111 


0.000 


0.058 


14 


44.99 


-267.836 


-507.013 


0.113 


-0.212 


11101 


0.033 


0.069 


15 


44.56 


-82.884 


-276.457 


-0.159 


-0.229 


01111 


0.031 


0.066 


16 


44.15 


-82.811 


-276.482 


-0.281 


-0.279 


11111 


0.031 


0.047 


17 


43.72 


-82.815 


-276.430 


-0.286 


-0.509 


11111 


0.021 


0.053 


18 


43.29 


-82.809 


-276.490 


-0.245 


0.037 


11100 


0.042 


0.019 


19 


42.87 


-155.561 


134.052 


0.562 


0.391 


11100 


0.004 


0.064 


20 


42.49 


-218.197 


104.703 


-0.117 


-0.039 


00111 


0.017 


0.006 


21 


42.07 


-218.168 


104.679 


-0.181 


0.022 


00111 


0.005 


0.002 


22 


41.65 


-218.141 


104.652 


-0.264 


0.103 


00111 


0.007 


0.012 


23 


41.62 


-271.528 


-48.440 


-0.097 


-0.117 


11111 


0.023 


0.032 


24 


41.20 


128.780 


-809.712 


0.454 


-0.599 


11111 


0.042 


0.078 


25 


40.76 


-215.215 


106.730 


0.131 


0.014 


11100 


0.020 


0.016 


26 


40.38 


-75.557 


73.994 


0.221 


0.100 


00111 


0.049 


0.011 


27 


40.23 


398.851 


-106.395 


0.553 


-0.164 


01111 


0.008 


0.082 


28 


39.96 


-51.367 


44.381 


0.102 


0.176 


00111 


0.017 


0.008 


29 


39.96 


-18.276 


-310.044 


-0.297 


-0.056 


00111 


0.028 


0.012 


30 


39.51 


-51.386 


44.441 


0.163 


0.116 


11111 


0.010 


0.022 


31 


39.51 


-18.438 


-310.014 


-0.127 


-0.103 


11111 


0.006 


0.015 


32 


39.09 


-18.481 


-310.026 


-0.018 


-0.104 


11111 


0.013 


0.010 


33 


38.66 


-18.555 


-309.907 


0.130 


-0.297 


01111 


0.010 


0.014 


34 


38.66 


-0.008 


0.022 


0.004 


-0.078 


11110 


0.003 


0.006 


35 


38.25 


0.000 


0.000 


0.000 


0.000 


11111 


0.000 


0.000 


36 


37.85 


-65.069 


-250.810 


-0.046 


-0.317 


00111 


0.010 


0.005 


37 


37.81 


-0.023 


0.029 


0.025 


-0.017 


11100 


0.019 


0.009 


38 


37.40 


-65.013 


-250.881 


-0.128 


-0.094 


11101 


0.001 


0.048 


39 


36.13 


-52.904 


-249.517 


0.822 


-0.185 


11110 


0.047 


0.024 


10 


35.72 


-52.870 


-249.514 


0.809 


-0.218 


11111 


0.054 


0.040 


11 


35.30 


-52.734 


-249.572 


0.576 


-0.064 


11111 


0.048 


0.032 


42 


34.88 


-52.602 


-249.609 


0.392 


-0.012 


11111 


0.034 


0.024 


43 


34.44 


-52.616 


-249.555 


0.443 


-0.058 


00111 


0.060 


0.023 



Column 1: spot ID number. 

Column 2: LSR velocity in km s — 1 (averaged over detected epochs). 
Column 3 & 4: best fit positions at the first epoch (in mas). 
Column 5 &; 6: best-fit proper motions (in mass). 
Column 7: detected epoch. '1' for detection, and '0' for non- 
detection. 

Column 8 & 9: RMS residuals from best fit proper motions (in 
mass). 



